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Using hard-loop (HL) effective theory for an anisotropic non-Abelian plasma, which even 
in the static hmit involves nonvanishing HL vertices, we calculate the imaginary part of 
the static next-to-leading-order gluon self energy in the hmit of a smah anisotropy and 
with external momentum parallel to the anisotropy direction. At leading order, the static 
propagator has space-like poles corresponding to plasma instabilities. On the basis of a 
calculation using bare vertices, it has been conjectured that, at next-to-leading order, the 
static gluon self energy acquires an imaginary part which regulates these space-hke poles. 
We find that the one-loop resummed expression taken over naively from the imaginary-time 
formalism does yield a nonvanishing imaginary part even after including all HL vertices. 
However, this result is not correct. Starting from the real-time formalism, which is required 
in a non-equilibrium situation, we construct a resummed retarded HL propagator with correct 
causality properties and show that the static limit of the retarded one-loop-resummed gluon 
self-energy is real. This result is also required for the time-ordered propagator to exist at 
next-to-leading order. 

PACS numbers: ll.lOWx, ll.lSBt, 12.38Mh 



I. INTRODUCTION 

Non-Abelian plasma instabilities d, 0| have received a lot of attention following the apparent 
failure of perturbative QCD at finite temperature to account for the rapid thermalization and strong 
collectivity deduced from the experimental results at the Relativistic Heavy Ion collider (RHIC) 
To leading order (LO) in the coupling and for small gauge field amplitudes, the dynamics of plasma 
instabilities are determined by the generalization of the hard-thermal-loop (HTL) gauge 



boson self-energy to anisotropic situations [Sl, l9|, llO|, lUfl ■ At LO, the scale associated with plasma 
instabilities is of the same parametric order as that of other collective phenomena such as the Debye 
mass, and consequently plasma instabilities strongly modify previous perturbative thermalization 



scenarios 



Much progress has been made in the study of the nonlinear evolution of non-Abelian plasma 
instabilities using real-time lattice simulations using hard-loop (HL) effective 
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theory [20!| (for classical field theory simulations see 

Also a number of analytical calculations have been performed. The leading order anisotropic 
gluon prop agator has been analysed in the temporal axial gauge in 0, l3, 11 1, and in covariant gauge 
in 25|, l26|] and 271]. Unlike the HTL gluon propagator, this anisotropic propagator contains non- 
integrable space-like poles, which signal the presence of instabilities in an anisotropic system. The 



HL gluon propagator has been used for perturbative studies of collisional energy loss [2 and the 
decay width of quarkonium bound states [30|, l3ll| in an anisotropic plasma, where these space-like 
poles do not present a fundamental problem. However, non-integrable singularities at zero frequency 
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occur in perturbative calculations of jet quenching and momentum broadening in the anisotropic 
quark- gluon plasma [sl, [H], where they appear to signal an enhancement of these observables 
compared to the equilibrium case. It was suggested in [3J| that these singularities are regulated by 
a non-zero imaginary part of the next-to-leading order (NLO) gluon polarization tensor. A partial 
one-loop result using the HL-resummed propagator but bare vertices was performed, suggesting 
that this was indeed the case. 

In this paper we shall show that the one-loop resummed expressions considered in [sl] do indeed 
yield a nonvanishing imaginary part after adding up all HL-resummed one-loop contributions in- 
cluding the nontrivial vertices present even at zero frequency. However, starting from the real-time 
formalism we show that these one-loop resummed expressions need to be evaluated using the static 
limit of the retarded propagator, but in the presence of plasma instabilities the retarded propagator 
is no longer obtained from the analytic HL propagator by a limiting procedure. Constructing and 
employing a resummed retarded HL propagator with correct causality properties we find that the 
NLO static gluon self-energy is real. 

In thermal equilibrium, the imagin ary part of the static gluon self energy would be required 



to vanish due to the KMS conditions [3j|. In [32| it was argued that, if the KMS conditions are 
violated in the anisotropic case, there might be a finite, discontinuous contribution. We note, 
however, that such a nonvanishing imaginary part to the self energy in the static limit would cause 
a fundamental problem with real-time perturbation theory, because the static limit of the time- 
ordered HL-propagator would become ill-defined. This can be seen as follows. In the limit of 
small anisotropics, we can use equilibrium distribution functions on soft HL-resummed lines, and 



therefore the time-ordered propagator is obtained from the retarded propagator using [35|] 

D'^{ko,k) = D'^^\ko,k) + 2n{ko)ilmD''^\ko,k), (1) 

where n{ko) is the Bose distribution function. The second term in this expression diverges for 
ko 0, if the imaginary part of the static self-energy is non-zero. 

The anisotropic NLO static gluon self-energy is obtained from the diagrams shown in Fig. [H In 



equilibrium, only static loop momenta have to be considered and no fermion loops contribute [36| 
at next-to-leading order, and we shall see that this is also true in the anisotropic case. Since the 
ghost self energy vanishes at leading order, the ghosts do not need to be resummed. The solid dots 
indicate leading order propagators and vertices which are obtained from the HL effective action 
(2q| . In contrast to the equilibrium situation, the HL vertices do not reduce to bare vertices in the 
static limit. The third diagram represents the HL counterterm which must be subtracted to avoid 
double counting. 





+ # t ^ 



FIG. 1: The diagrams that contribute to static gluon self energy. All lines correspond to gluon propagators. 
The dots indicate hard loop propagators and vertices. The cross denotes the counterterm. 



At zero temperature, field theory can be formulated covariantly. At finite temperature, covari- 
ance is broken by the vector = (1,0,0,0) which specifies the rest frame of the thermal system. 
For anisotropic systems we need (in the simplest case) one additional vector to specify the direction 
of the anisotropy. We consider the case in which there is one preferred spatial direction along which 
the system is anisotropic (in planes transverse to this vector the system is isotropic). In the context 
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of heavy ion collisions, we can take this direction to be the beam axis (z) along which the initial 
expansion occurs. 

In [27|| we have given the complete analytic result for the integrand corresponding to the diagrams 
in Fig. [H Although this result is relatively compact, the evaluation of the integrals would be a 
formidable task. In order to simplify the calculation, we choose a particular orientation for the 
external momentum. We calculate the part of the NLO static gluon self energy related to Weibel 
instabilities for an anisotropic plasma, in the limit of small anisotropy, with the external momentum 
parallel to the direction of the anisotropy. 

This paper is organized as follows. In section [U we define our notation. In section IIIII we 
present the calculation of the integrand. In section [IV] we discuss the analytic structure of the 
HL propagator. In section |V] we calculate the NLO contribution using the analytic propagator 
obtained by resumming the HL gluon self-energy. In section |Vl] we construct a retarded HL- 
resummed propagator and identify an extra contribution that cancels the imaginary part obtained 
in section [Vl In section [VTTl we present our conclusions. 



II. NOTATION AND INGREDIENTS 

Throughout this paper we will frequently use the indices {k, q, r} to denote momentum argu- 
ments. We also use Latin letters /,•••} to denote spatial indices, with the exception that the 
indices k, q, r are reserved and used exclusively to denote spatial momenta. Four-momenta will 
be denoted by capital letters. The external momentum in Fig. [T]will be called q, and the internal 
momenta are k and r = —k — q. We choose q = (0, 0, q). 



A. Distribution Functions 

For a non-Abelian plasma with Nc gluons and A'^^ quark flavors, we define an effective isotropic 
distribution function by 

/iso(p) = 2Nf[nip) + hip)] + 4NcnS{p). (2) 
In thermal equilibrium, we have 

We define the Debye mass from the equilibrium distribution^ 



Following pj|, we can construct an anisotropic distribution from any isotropic distribution of 
the form /iso(p^) by writing 



where > —1 is the anisotropy parameter. A value ^ > corresponds to a contraction of the 
distribution and > ^ > — 1 corresponds to a stretching of the distribution. In the following we 



[1] In Ref. [01 Eq. (g} differs by a factor of 2 and Eq. ^ differs by a factor 1/2. The definition of the Debye mass 
is the same. 
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shall consider only deformations of the equilibrium distribution. For nonzero ^, the parameters T 
and ^ of course lose the usual meaning of temperature and chemical potential. In this paper we 
restrict ourselves to weakly anisotropic systems for which |^| <C 1 and shall calculate only to leading 
order in ^. 



B. Vertices 

Momenta are taken to be incoming. We give only the tadpole form of the 4-point vertex, since 
that is the only 4-point vertex we will need. The bare vertices are: 

(ro)::,^,' = ^ff/a6c^r' (6) 

(MoTak^^iQ, -Q, K, -K) = 2g^ CA^ab iVC"'^ 
The notation for the HL vertices is given in 2^, 21\- We define: 



1 



,p (2vr)2 
The 2-point function is: 



po=p 



; P^:=il,f); I, ■■= [ ^ 



p 



P-Q 



The 3-point function is: 

■■= ^9fabcr^'' (9) 

\P-kP-q P-qP-r 

We need only the tadpole form of the 4-point vertex which has the form: 

<6cf(Q. -Q^ -K) := 2g^ CA^ab M^^^^^^Q, -Q, K, -K) (10) 



M'"'^''{Q, -Q, K, -K) := -llfsP'^P" P^P° 



K^P ■ Q - Q^P ■ K 



P-K P-Q{{P-Kf-{P-QY) 



These vertices satisfy the Ward identities: 

i^^r^^^iT, Q, R) = w^{Q) - n^^(i?) (11) 

KxM'^'^^Q, -Q, K, -K) = -2T'''"'{K, Q, -K - Q) 
KxK„M^'''^''{Q, -Q, K, -K) = 2U^''{-K - Q) - 2I^^''' {Q) 



C. Leading order hard- loop self energies and propagators 

In Feynman gauge, which we use, the relation between the polarization tensor 11 and full (bare) 
propagators D (D^) is given by the relation 

d^Xk) = {D%r' - n^. = -{g^uK^ + n^.)- (12) 
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The HL gluon self-energy is gauge invariant and satisfies the usual Ward identity: K^^H^y = 0. 
Consequently, we only need to calculate the spatial components. 

In the isotropic (HTL) case, the spatial self energy has two independent components which are 
called the transverse and longitudinal parts. Eq. (fT2]l defines an analytic propagator from which 
the retarded (advanced) propagator is obtained as 

D^-^/-d-[K) = D^,{k,±ie,k). (13) 

For anisotropic systems the self energy can be generally decomposed into four independent struc- 
ture functions. For the special case of n^i,(go,9) with q = (0,0, g), there are only two independent 
structure functions. We define the general projection operators using the vector [37|] 

4 := n%k) = - k'k^/k^)5^^ (14) 

which satisfies fc* = 0. Using this vector we construct the projection operators for Ilij{K): 

plk _ c _ kjkj 2k _ kj^j p3k _ _kikj_ _ , e 

~ A;2 ' ~ A;2 ' ~ k'^ ^1 ' 

P^J' = k,n^ + hn';, P^ = P^-P^, (15) 

and similarly for Ilij{R). The projection operators for JIij{Q) are defined as: 

P^ = 6^,-6si6s,; P!f = 63,63j (16) 

The orthogonality relations satisfied by these projectors are given in Appendix A. Using these 
definitions, the self energies can be decomposed as: 

nf,=^^>. + ^^^, (17) 



nf^. = Plj^ak + PiJ" h + p!- Ik + Ptj" h 



The complete LO results for the HL self energy components a^, (3k, 7fe, 6k have been given in [ul] 
for the anisotropic distribution function jS]). 



D. Static propagator 

The static propagator has (space-like) poles at go = and values of q which mark the boundaries 
between stable and unstable modes. These are determined by the static parts of au and 7^, which 
are given by: 

(^k = -^{^-nl)m\)- ^k = ^nlm% (18) 



For the (Feynman gauge) HL propagator Dij{q) with q = (0,0, g) we obtain, in the static limit, 

1 pi, 1 

q2 rj g2 _|_ Q,^ 



DM = -tf-2-Plf:;T^- (19) 



For ^ > 0, the second term has a space-like pole at = _q,^ > g, which corresponds to the 
magnetic Weibel instability [s^]. More precisely, all modes for which < q^ < {—aq) exhibit 
exponential growth with a growth rate 7(g) that approaches zero as g^ approaches or —Og (see 



Fig. [21). The result for the general static propagator D^^{k) is given in Ref. [27|. To leading order 
in ^, it has space-like poles at fc^ = —q-^ and fc^ = — (a/c -|- 7^) determined by 
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III. INTEGRAND FOR a, AT NEXT-TO-LEADING ORDER 



To obtain Imanio we use (fTTll and write 



Im anlo = ^Pl^iq) Im U,,{q) (20) 

In order to calculate lmW^{q) we need the integral corresponding to the first two diagrams in 
Fig. [H Since the imaginary-time formalism is not available in the nonequilibrium case we are 
considering, our starting point is the (resummed) one-loop expression for the retarded self-energy 
as given by the real-time formalism. We use the fact that the vertices are pure real in the static 
limit (see Appendix B) and obtain: 



W^iq) = -g'CA 



A" J (27r) 



(21) 



where we have defined: 



= ((ro + r)^-prA'(r° + r)^'^--'p;!^r: 
Hli = ((M° + uyi^'^p^^] 



(22) 

go=fco=0 



go=fco=0 

D^y^{K) = (1 + 2nl^{ko)){Dl'\K) - Df^\K)) 

DT'iR) = (1 + 2nl^{ro)){Dl^\R) - Df"{R)) 

The order of the momentum variables is [K, Q, R) for the 3-point functions and (Q, —Q, K, —K) for 
the 4-point functions. The repeated indices n and m indicate sums over the projection operators 
in (fTSl) . and the corresponding mode functions in the propagator. These functions are given in 



27|| and we do not reproduce them here because, as we shall see below, not all components of the 
propagator contribute at leading order. 

To lowest order in ^ we can use the equilibrium distribution function in the third and fourth lines 
in (f22]l . The dominant contribution to the integral can be extracted by using (l-|-2nfq(A;o)) 2T/kQ 
which gives: 

^ 1 9 ^ „ /' d^k f dkn 1 , , 

n-(,o = 0,.1 =-,9^C,tJ^^ j^- (23) 

■[Fii^-i {Dl^\K)Dlf{K + Q)- Dt'{K)D^:i^{K + Q)) 

+/f^.i(D-*(K)-Z)rW)] 1,0=0 

Since D^^^ = (D"-'^'")* , this expression appears to be manifestly real. However, it is not well- 
defined, because in the anisotropic case, there are space-like poles in the propagators at ko = 0. 
We discuss this issue in the next section. 



IV. ANALYTIC STRUCTURE OF THE HL PROPAGATOR 



The a naly tic structure of the anisotropic HL propagator defined by (fT2]l has been discussed in 
Refs. [10|, lull- The poles in the propagator at real frequencies and real wave vectors determine the 
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dispersion laws of propagating modes. Poles at real frequencies and imaginary wave vector describe 
(frequency- dependent) screening effects. In the anisotropic case, there appear also poles at real 
wave vector and imaginary frequency corresponding to plasma instabilities. As an example, Fig. [2] 
displays the location of all these poles in the propagator Da{(jJ, 0, 0, q) in a plot of cu^ versus /c^ for 
3 values of ^. 




FIG. 2: Poles in the propagator component Da{uj,q) for q parallel to the direction of anisotropy, for 
anisotropy parameter ^ — —0.5 (prolate momentum distribution), ^ = (isotropic case), and ^ = +5 
(oblate case, with space-like pole at zero frequency). 



The analytic HL propagator Dfj_^{ko, k) with complex ko but real k generally has time- like poles 
on the real axis at ko = uj{k) > /c, a logarithmic branch cut between ko = ±k, and, for sufficiently 
small k (with details depending on the direction of k and the polarization) poles at ko = ±i\'y{k)\, 
where \^{k)\ is the growth rate of unstable modes (not to be confused with the structure function 7^ 
in the decomposition ofIlij{K)). When these additional poles are present, using Cauchy's theorem 
and an integration contour in complex frequency as depicted in Fig. [3j the spectral representation 
of the analytic propagator has the form: 

^) = r ^fco Djko + ie, k) - Djko - ie, k) _ sp ResZ)(±i|7(fc)|, fc) 
y_oo27ri ko-uj ^ ±z|7(A?|)-w 



The last term in ([24ll is produced by the poles at ±z7(fc). In the isotropic case this term is absent 
and we recover the usual spectral representation. 
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FIG. 3: Contour used to derive Eq. i|24p and Eq. (|25| . Crosses and thick lines indicate poles and cuts, 
respectively, in the anisotropic HL propagator. The poles on the imaginary axis appear at sufficiently small 
k and disappear from the physical sheet at higher k. In the case of Eq. l(25|) . there may be an additional 
pair of poles on the imaginary axis from the second propagator in the one-loop diagram that needs to be 
encircled clockwise. 



However, whenever the analytic propagator D{z, k) has poles on the imaginary axis correspond- 
ing to plasma instabilities, D{uj + ie, k) does no longer yield the retarded propagator, since it does 
not satisfy the condition that J duje~^^^ D{uj + ie,k) vanishes for t < 0. This point will be discussed 
in section |Vl] where we present a modification of the formalism in which the retarded propagator has 
the correct causal structure. In section|V]below, we consider the consequences of simply ignoring the 
issue, and calculating the next-to-leading order polarization tensor by defining retarded/ advanced 
propagators as usual through (fT3ll . 



V. NEXT-TO-LEADING ORDER CALCULATION OF lma(g=0,g) USING 
UNMODIFIED HL RESUMMATION 



We begin by following the strategy that would work in standard HTL perturbation theory: We 
can rewrite the /cQ-integral in ([231) using a sum rule technique (see, for example, ^^''^ ^ 

away from the real axis, we can use Cauchy's theorem to write: 

dkf) X{ko + ie, k) — X{kQ — ie, k) sr-^ ResX(i7s(fe), k) 
2'Ki ko — k^ ^ i"fs{k) — k'^ 

dko X'-^^ko, k) - X'"^''{ko, k) ResX{i-fsik), k) 

27ri /co — feg i^sik) — k'Q 



oo 



where the contour C is analogous to the one in Fig. [U except that there will be generally two pairs 
of simple poles appearing symmetrically on the imaginary axis, from D{ko,k) and D{kQ,r}, which 
are denoted collectively by 7^. In the last line we have tentatively assumed the validity of Eq. (fT3]l . 
which will be reassessed in the next section. 

We have to consider X{ko, k) = hH Dn{K) and X{ko, k) = 1^ Fi^mDn{K)Dm{K + Q). Taking 
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the limit /cq ^ 0"'" = lim£_>o lim^^o V + i^, Eq. ([25!l gives: 
Im 



°° d/to 1 



— oo 
oo 



Z/;^ . (Z^r (^) - = -Hli Im lim Z)-* (ET) 

Im / ^ ^ ^ F^i, z {Dl^\K)DZ\K + Q) - Dl'^{K)D^^\K + Q)) 

-F^i, Im lim Z?-*(i^)D-*(i^ + Q) 

fcO-*0+ 



(26) 



The second term on the right hand side in (|25|1 does not contribute to the imaginary part 
since the residues ResX(z7(A;), /c) are purely imaginary and odd in 7. This means that 
Y^'ResX{i^s{k),k) / {i'^s{k) — k'o) is purely real when fcg approaches the real axis, since for each 
value of 7 there is another with reversed sign. 
Substituting ([261) into (f23l) we obtain: 



(27) 



(27r)3 



[ F^i, Im , lim , D':\K)DZ\K + Q) + Im , lim , D^:\K) 



fcO-+0+ 



go=o 



The integral in (f27l) could have an imaginary part coming from the space-like poles of the static 
internal propagators at fco = when k is just at the boundary of the instability domain. This 
observation led to the conjecture in Ref. [32l| that this imaginary part is in fact nonzero, and that it 
regulates the singularities in the momentum broadening coefficient that are produced by the space- 
like poles in the static HL propagator. This conjecture was supported by a partial calculation of 
the integral in ([27l) . in which only the tadpole diagram with a bare vertex was included. Below we 
calculate the complete result including all HL vertices which are nontrivial in the anisotropic case, 
in contrast to their HTL counterparts. 

To leading order in ^, the poles in the fc-propagator occur at /c^ = —a^ and A;^ = — (a^ +7a;), or 
A; ~ + 0(6- Therefore, we need only the part of this propagator that is leading order in ^ 

for momenta k ~ mo- It is easy to show that only the spatial components contribute. Since we 
are working on the LO mass shell q = \/^/^m£), the same conclusion holds for the r-propagator. 
Assuming ([H]), we obtain: 



lim Dlf{K) 



LO 



(28) 



lim 



F 



2k 



1 



I 't '} 



3k 



1 



kl-k"^ + i sgn(/co)e ''^ kl - {k'^ + au) + i sgn(/co)e 



.5fe 



kl - (fe2 + ak + -fk) + i sgn(A;o)e 
The real and imaginary parts are: 

Re [ , lim , Dlf{K) ] ^ -Pf^^^ - Pf^^h{k) - Pf^^ h{k) 



lim 

fco^O+ 



(29) 



^Im [ , hm^ D^{K) ] ^ -.vr[/^f ^(/g-i) + Pfl'dif,')] 



h{k) 



fc2 + ak 



; h{k) 



k"^ + ak + 7fc 
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Insertion into (I27|l gives: 
1 



n,mG{3,5} 
nG{3,5} 



(30) 



We obtain Imanio by substituting (|30]) into (|20|) . We give below the result for the integrand for 
the tadpole graph. The corresponding expression for the bubble graph is straightforward to obtain, 
but considerably longer. 



Ima 



nlo 



tp 



(2^ 



(31) 



[ - M^'" 5{f^,^) - 4(fc2 + 3^3^)^ 5{f-,^) - 4(A;2 - ?,q{k; + nlq/2))^5{f~^^)\ 
In order to perform the integrations, we choose the following representation: 

n = (0,0,1) (32) 

g = (0,0,g) 

k = {k sin 0, 0, k cos 9) 

pi = (sin 6*1 cos (pi , sin 6i sin 0i , cos 9i ) 

P2 = (sin 02 cos (/>2, sin 02 sin 02, cos 62) 

The vector n gives the direction of the anisotropy. The vectors pi and p2 are used to calculate 
the HL vertex components M*'-'-^ and P-^'P-''. We will write cosO = x, cos^i = xi, cos02 = X2. 
We work on the LO mass shell q = m/j. Using the fact that each term in (f3T]l contains a 

delta function which gives k ~ \/^mD, we scale /c to obtain a dimensionless variable k which is 
defined by: k = {\/C/3mD) k. Using this coordinate system, all angular integrals in the HL vertices 
can be calculated analytically. Next, we use the delta functions to do the polar integral over the 
variable x = cos 9. This produces constraints on the k integration. Finally, the resulting /c-integral 
can be done by numerically. The full integrand produces a huge number of terms, many of which 
are divergent. When all terms are combined, all divergences cancel, which provides a check of our 
algebra. Some details of the calculation are given in Appendix [Bl The final result is: 



Im(aq)nio = 9^ NcVnoT 



(33) 



I • 3.77... 

This finite result confirms the conjecture of Ref. [3^ that a complete one-loop evaluation of Eq. (f27ll 
using full HL vertices and HL propagators would lead to a nonvanishing imaginary part to the static 
gluon polarization tensor that may be used as a regulator for the space-like poles in the static gluon 
propagator of an anisotropic plasma. 

However, as already mentioned in the previous section, the finite result in ([33]) is produced by 
using HL retarded and advanced propagators which do not have the correct causal structure. In 
fact, it is easy to see that there is a fundamental problem with the above calculation. The sign of 
the imaginary part in ([331) is determined by having taken the static limit /cq = + — > 0+ in ([271l . 
Taking the static limit e.g. instead by /cq = —fl + would produce the opposite sign. In the next 
section we discuss how to define retarded and advanced propagators that have the correct causal 
structure. We then show that using these propagators produces a pure real result. 
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FIG. 4: Complex frequency contour required to obtain the retarded propagator from the analytic propagator. 



VI. RETARDED HL PROPAGATOR AND CORRECTED HL RESUMMATION 

Using the analytic propagator D obtained from Dyson summation of the HL polarization tensor 
according to Eq. p^ . we can construct a retarded propagator 

D'^^t{t,k)= [ duje-''^^D{uj,k) (34) 



I a 

with the property that D'^'^*{t < 0, A;) = by choosing a complex contour Cr that runs along the 
real axis and circumvents any poles at positive imaginary oj by going over them, as shown in Fig. [H 
Moreover, for t > this propagator correctly describes the exponential growth of the unstable 
modes, which an integration over just the real axis would have missed. Restricting uj to real values 
again, we can take the effect of the contour Cr into account by defining 

Z)-*(a;, k) = Dr^iu; + ie, k) - ^) e(^„(g) _ e) (35) 

UJ + ie- ijn{K) 

where Dn represents a component of the analytic propagator. For k large enough so that all modes 
are stable, all 7„(A;) are negative and this definition of the retarded propagator coincides with the 
usual one. But when k becomes small enough for unstable modes to appear, the corresponding 
poles are effectively subtracted. 

For example, in the small-^ limit we have 

a{0,k) = -(C/3)m|, 008^6* + 0(^2) with cosO = k^/k 

Ah 

7(fc) = V(^'-(^/3)^l)Cos2 + O(a) 

-> Aik 

ResDaii-f,k) = ^ + O(0 (36) 

so that the static limit of the retarded static a-propagator responsible for the Weibel instabilities 
now reads 

Dc,iO,k)=— a"*" — 2 a I r,fc2\ for A; > yO^mu cos 6* 

k'^ — (4/3)m2) cos^ d + 



Da{0,k)= 0(^°) for k < ^/^mn cos 0. (37) 

In particular, the imaginary part proportional to 6{k'^ + a{0,k)) is cancelled (irrespective of the 
sign of the infinitesimal rj in uj = rj + ie) , and the same is true for all space-like poles in D (and 
also does not require the simplifying assumption of momentum parallel to the anisotropy direction 
used in our explicit calculations). Thus the entire contribution to the imaginary part calculated in 
section |V] disappears. 



12 



VII. CONCLUSION 

Starting from the one-loop expressions of the real-time formalism we have found that these 
require retarded and advanced propagators, which in the anisotropic case after resummation of 
the HL gluon self-energy are no longer given as boundary values of the analytic propagator in the 
limit of real frequencies. Using the modified prescription (|35p we found that the imaginary part 
of the static gluon self-energy vanishes also at HL-resummed one-loop order. This disproves the 



conjecture of Ref. [3j| that such an imaginary part would be generated, and then could be used to 



regulate the non-integrable singularities encountered in p erturbative calculations of jet quenching 



and momentum broadening in an anisotropic plasma [32, 13; 



In Ref. [33|| it has been pointed out that the space-like poles present a problem only together with 
the singularity of the Bose- Einstein distribution function, and that, as an alternative regularization, 
it would be natural to assume a lower cutoff on the frequencies because of the inherently nonequi- 
librium nature of the problem, involving characteristic time scales ^ (q£,T)~ ^. This resolution to 
the problem of non-integrable singularities in the calculations of Ref. U 33 1 in fact would lead to 



a larger enhancement of the effects of anisotropy on jet quenching and momentum broadening than 
those resulting from a nonzero imaginary part of the gluon self-energy of order (/^Tm^i^/f. Our 
results evidently lend support to this possibility. 
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APPENDIX A: PROJECTION OPERATORS 

The orthogonality relations satisfied by the operators given in section Hi CI are: 

PilPij = ; PiiP^j = ; Pii^tj = ; Pu^tj + Pu^ij = -^4' ^^^^ 

T}2 p2 p2 . p2 p3 n ■ p2 p4 , p4 p2 p4 

p3 p3 p3 . p3 p4 n 

p2 p5 n • p3 p5 n • p4 p5 p4 . p5 p5 p5 

h — k P^ — h k P^ — h P^ — n 

Tr = 2 ; Tr = Tr = Xr = 1 ; TrP^ = 

TrP^'i = 2 ; TrP^? = 1 ; pi? • P^? = Q (A2) 

APPENDIX B: SOME DETAILS OF THE CALCULATION IN SECTION [V] 

1. HL vertex components 



The only vertex components that we need are M^^i^ and r^^'F'J'. Using ©, ([lO]) and ([32|), these 
vertex components can be rewritten as: 
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1 



27r rl 



M''" = - I d4>, I dxi ■ Im, (B1) 



1 



6x 3kx 
Im = ^ , ■ , 1- y 



Xxl + Vl — COS — x\xi kxx\ — x\ + k\/l — x'^ COS — x\xi 

3kx 3 



+ ^ 



kxxf + xf + k\/l — 2;2 cos (pi ^J\^^x\x\ kxx\ — x\^ k\J\ — x"^ cos -y/l — x^xi 
3 



fcxxf + + fcVl — cos 01 y^l — x^xi 
riiip^-« = 1 ^ ^ y dxi y dx2 ■ /rr, (B2) 



3X^ ( Y^l - xf Y^l - X2 cos ((/>l - 02) + XlX2 

Irr = 



\/l — x^ Y^l — x^ cos (pi + xxij ^Vl — x^ Y^l — X2 cos 02 + a;x2 

3x(fcx + 1) (^1/1 - xf 1/1 - X2 cos (01 - 02) + 2;iX2^ ^ 
k\/\ — X^y^l — x\ COS 01 + A:XXl + Xi^ (^\/ 1 — x'^ -s/ 1 — x'2 cos 02 + XX2 
3x(fcx + 1) ( Y^l — Xi Y^l — X2 COS (01 — 02) + X1X2 ' ^ 



\/l — x^Y^l — x^ COS 01 + xxi J ( A;-v/l — x^ y^l — x^ cos 02 + A;xx2 + X2 



3(A;x + 1)^ ( v^l - xf y^l - X2 COS (0i - 02) + X1X2 



k\J\ — x'^y^l — x'f COS 01 + A:xxi + 2;i ) ( k\^l — x^ y^l — Xg cos 02 + fcxx2 + X2 



Using the identities cos(0i — 02) = cos 0i cos 02+sin0i sin 02 and sin^ 0i = 1— cos^ 0i the azimuthal 
integrals in Eqs. (IBip and (IB2p can be rewritten in the form (with the same expressions for the 02 
integrals) : 



1 _ sgn(^)e - B^) ie (52 - A^) 



2tt A + B cos 01 - 52 V52 _ ^2 

d0i cos 01 _ Asgn(.4)9 (^2 - ^2^ iAQ {B^ - A^) l_ 

2tt A + B cos 01 ~ BVA^ - B^ ^ BVW^^^ ^ B 

d(l)i cos2 01 _ sgn(^)9 (^2 - B^) A^ _ iQ {B^ - A^) A^ _ A_ 

2tt A + B cos 01 ~ BWA^ - B^ B^VW^W ~ 

dcPi cos3 0i _ sgn(^)e (^2 - 52) ^3 iQ{B^-A'^)A^ 2A^ + B'^ 



(B3) 



2ti A + B cos 01 SV^2 _ 52 5V52 _ ^2 2^3 

After performing the azimuthal integrals, we are left with the polar integrals over the variables 
xi and X2 with limits determined from the theta functions in Eqs. (IB3p . In every case, the imaginary 
part contains the integral of an odd function over a symmetric interval and therefore, as is the case 
for the self energies, the imaginary part of all components is zero in the static limit. This result 
corresponds to the fact that, in the static limit, the bare propagator has no poles. The xi- and X2- 
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integrals can be done analytically, yielding 



m2 (x2 - 1) (8x2 + ^ /g^s _^ 22x) /c + Sx^ + 7 U 

= 0, r^^r^^ = ^ (B4) 

3(P + 2xk + 



2. Extraction of the imaginary part 



We substitute the results (|B4I) into the integrand and do the x-integrals. In order to extract the 
imaginary part, we write (with corresponding expressions for fsir) and f^{r)): 

Using Eqs. (fTSl) we can rewrite these expressions: 

1 3 

(B6) 



k'^ + ak-ie ^m\,{x - ie - x^^){x + ie - x'^^) 

1 _ P + 2xfc + 1 

r"^ + ar - ie k'^^m\){x + ie - x'^^){x - ie - xj^) 
1 3 



fc2 + Ofc + 7fc - ie 2^mjj{x - ie - Xg^)(x + ie - x+^) 
I 3 (P + 2xk + 1^ 



with 



+ ar + jr - ie 2A;2^m|)(x + ie - x'^j.){x - ie - Xg^) 
^3k = -k; x'^^ = k (B7) 



X 



3^ — — K ; xjj. = — when /e < i 

+ r - ~k'^ + 2 ~ 
= —A; ; Xg^ = = — when fc > 1 



3/c 



\/PTl . \/A;2 + 1 



V2 ' ^'^ V2 



X 



2P + v^2A:4 _ 2A;2 + 1 + 1 -2p + \/2A;4 - 2A;2 + i _ i 

2^ ' 2i 



Substituting (IBSP and (IB6j) into ([3T]) the x-integrals can be done analytically. After a partial 
fraction expansion, the integration of each term produces either a logarithm, or something that 
is pure real. The imaginary part of each logarithm gives rise to a theta function, restricting the 
remaining fe-integrals to intervals 0</c<l,l<A;<2, orO<A;<2. We note that particular care 
must be taken because of the double pole at x^^ = Xg^ = —k that occurs when k < 1. We regulate 
this term in the usual way by writing: 

I d I , , 

(B8) 



i^-Hk){x-Hk) dMx-x^^-M 



and taking M to zero at the end of the calculation. 
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Once all the x-integrals are done, the remaining integral over k can be done numerically. There 
are many terms that are divergent at k = 0, k = 1, and/or k = 2. Summing all contributions yields 
the result given in Eq. (f33]) . 
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